PaCS-Toolkit: Optimized Software Utilities for Parallel Cascade Selection Molecular Dynamics (PaCS-MD) Simulations and Subsequent Analyses.

Parallel cascade selection molecular dynamics (PaCS-MD) is an enhanced conformational sampling method conducted as a "repetition of time leaps in parallel worlds", comprising cycles of multiple molecular dynamics (MD) simulations performed in parallel and selection of the initial structures of MDs for the next cycle. We developed PaCS-Toolkit, an optimized software utility enabling the use of different MD software and trajectory analysis tools to facilitate the execution of the PaCS-MD simulation and analyze the obtained trajectories, including the preparation for the subsequent construction of the Markov state model. PaCS-Toolkit is coded with Python, is compatible with various computing environments, and allows for easy customization by editing the configuration file and specifying the MD software and analysis tools to be used. We present the software design of PaCS-Toolkit and demonstrate applications of PaCS-MD variations: original targeted PaCS-MD to peptide folding; rmsdPaCS-MD to protein domain motion; and dissociation PaCS-MD to ligand dissociation from adenosine A2A receptor.

Delineating the conformational landscape of the adenosine A2A receptor during G protein coupling.

G-protein-coupled receptors (GPCRs) represent a ubiquitous membrane protein family and are important drug targets. Their diverse signaling pathways are driven by complex pharmacology arising from a conformational ensemble rarely captured by structural methods. Here, fluorine nuclear magnetic resonance spectroscopy (19F NMR) is used to delineate key functional states of the adenosine A2A receptor (A2AR) complexed with heterotrimeric G protein (Gαsß1γ2) in a phospholipid membrane milieu. Analysis of A2AR spectra as a function of ligand, G protein, and nucleotide identifies an ensemble represented by inactive states, a G-protein-bound activation intermediate, and distinct nucleotide-free states associated with either partial- or full-agonist-driven activation. The Gßγ subunit is found to be critical in facilitating ligand-dependent allosteric transmission, as shown by 19F NMR, biochemical, and computational studies. The results provide a mechanistic basis for understanding basal signaling, efficacy, precoupling, and allostery in GPCRs.

Kinetic Selection and Relaxation of the Intrinsically Disordered Region of a Protein upon Binding.

Here, we investigate the association and dissociation mechanisms of a typical intrinsically disordered region (IDR), transcriptional activation subdomain of tumor suppressor protein p53 (TAD-p53), with murine double-minute clone 2 protein (MDM2). Using a combination of cycles of association and dissociation parallel cascade molecular dynamics, multiple standard molecular dynamics (MD), and the Markov state model, we were successful in obtaining the lowest free energy structure of the MDM2/TAD-p53 complex as the structure closest to the crystal structure without prior knowledge of the crystal structure. This method also reproduced the experimentally measured standard binding free energy, and the association and dissociation rate constants, requiring only an accumulated MD simulation cost of 11.675 µs even though that actual dissociation occurs on the order of seconds. We identified few complex intermediates with similar free energies; yet TAD-p53 first binds MDM2 as the second lowest free energy intermediate kinetically with >90% of the flux, adopting a conformation similar to that of one of these few intermediates in its monomeric state. Even though the mechanism of the first step has a conformational-selection-type aspect, the second step shows induced-fit-like features and occurs as concomitant dehydration of the interface, side-chain π-π stacking, and main-chain hydrogen-bond formation to complete binding as an α-helix. In addition, dehydration is a key process for the final relaxation process around the complex interface. These results demonstrate that TAD-p53 kinetically selects its initial binding form and then relaxes to complete the binding.

Dissociation Process of a MDM2/p53 Complex Investigated by Parallel Cascade Selection Molecular Dynamics and the Markov State Model.

Recently, we efficiently generated dissociation pathways of a protein-ligand complex without applying force bias with parallel cascade selection molecular dynamics (PaCS-MD) and showed that PaCS-MD in combination with the Markov state model (MSM) yielded a binding free energy comparable to experimental values. In this work, we applied the same procedure to a complex of MDM2 protein and the transactivation domain of p53 protein (TAD-p53), the latter of which is known to be very flexible in the unbound state. Using 30 independent MD simulations in PaCS-MD, we successfully generated 25 dissociation pathways of the complex, which showed complete or partial unfolding of the helical region of TAD-p53 during the dissociation process within an average simulation time of 154.8 ± 46.4 ns. The standard binding free energy obtained in combination with one-dimensional-, three-dimensional (3D)- or Cα-MSM was in good agreement with those determined experimentally. Using 3D-MSM based on the center of mass position of TAD-p53 relative to MDM2, the dissociation rate constant was calculated, which was comparable to those measured experimentally. Cα-MSM based on all Cα coordinates of TAD-p53 reproduced the experimentally measured standard binding free energy, and dissociation and association rate constants. We conclude that the combination of PaCS-MD and MSM offers an efficient computational procedure to calculate binding free energies and kinetic rates.

Vibrational Energy Transfer from Heme through Atomic Contacts in Proteins.

A pathway of vibrational energy flow in myoglobin was studied by time-resolved anti-Stokes ultraviolet resonance Raman spectroscopy combined with site-directed mutagenesis. Our previous study suggested that atomic contacts in proteins provide the dominant pathway for energy transfer while covalent bonds do not. In the present study, we directly examined the contributions of covalent bonds and atomic contacts to the pathway of vibrational energy flow by comparing the anti-Stokes resonance Raman spectra of two myoglobin mutants: one lacked a covalent bond between heme and the polypeptide chain, and the other retained the intact bond. The two mutants showed no significant difference in temporal changes in the anti-Stokes Raman intensities of the tryptophan bands, implying that the dominant channel of vibrational energy transfer is not through the covalent bond but rather through van der Waals atomic contacts between heme and the protein moiety. The obtained insights contribute to our general understanding of energy transfer in the condensed phase.

Ligand binding to anti-cancer target CD44 investigated by molecular simulations.

CD44 is a cell-surface glycoprotein and receptor for hyaluronan, one of the major components of the tumor extracellular matrix. There is evidence that the interaction between CD44 and hyaluronan promotes breast cancer metastasis. Recently, the molecule F-19848A was shown to inhibit hyaluronan binding to receptor CD44 in a cell-based assay. In this study, we investigated the mechanism and energetics of F-19848A binding to CD44 using molecular simulation. Using the molecular mechanics/Poisson Boltzmann surface area (MM-PBSA) method, we obtained the binding free energy and inhibition constant of the complex. The van der Waals (vdW) interaction and the extended portion of F-19848A play key roles in the binding affinity. We screened natural products from a traditional Chinese medicine database to search for CD44 inhibitors. From combining pharmaceutical requirements with docking and molecular dynamics simulations, we found ten compounds that are potentially better or equal to the F-19848A ligand at binding to CD44 receptor. Therefore, we have identified new candidates of CD44 inhibitors, based on molecular simulation, which may be effective small molecules for the therapy of breast cancer.